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Abstract — The link scheduling in wireless multi-hop networks 
is addressed. Different from most of work that adopt the protocol 
interference model which merely take consideration of packet 
collisions, our proposed algorithms use the physical interference 
model to reflect the aggregated signal to interference and noise 
ratio (SINR), which is a more accurate abstraction of the real 
scenario. We first propose a centralized scheduling method based 
on the Integer Linear Programming (ILP) and resolve it by an 
approximate solution based on the randomized rounding method. 
The probability bound of getting a guaranteed approximate 
factor is given. We then extend the centralized algorithm to a 
distributed solution, which is favorable in wireless networks. It 
is proven that with the distributed scheduling method, all links 
can transmit without interference, and the approximate ratio of 
the algorithm is also given. 

I. Introduction 

In wireless multi-hop networks, concurrent transmissions 
that share a common channel may cause interference, and 
if too many devices ttansmit simultaneously, the interference 
will prevent an intended receiver from receiving the signal, and 
causes message loss. On the other hand, if too few nodes trans- 
mit at the same time, valuable bandwidth is wasted and the 
overall throughput may degenerate. Hence, the classic problem 
faced by the MAC layer scheduling protocols is to select an 
appropriate set of devices for concurrent ttansmissions, so 
that the interference does not cause message loss. In slotted 
wireless multi-hop networks, a natural and important goal of 
scheduling algorithms is to maximize the total throughput with 
the interference restriction. 

Hence, an accurate model of interference is fundamental in 
order to derive theoretical or simulation-based results. There 
are two main interference models used, namely the protocol 
interference model and the physical interference models [1]. In 
the protocol model, a communication from node u to node v is 
successful if no other node within a certain interference range 
from v is simultaneously ttansmitting. Due to its simplicity 
and to the fact that it can be used to mimic the behavior of 
CSMA/CA networks such as IEEE 802.11, this model has 
been widely used in the literature. However, it doesn't reflect 
the advanced physical layer technologies that allow concurrent 



multiple signal reception. In the physical interference model, 
a communication between nodes u and v is successful if the 
Signal to Interference and Noise Ratio (SINR) at receiver v is 
above a certain threshold. This model is less restrictive than the 
protocol interference model, and thus higher network capacity 
can be achieved by applying the physical interference model. 

Clearly, the interference in the protocol model is a tremen- 
dous simplification of the physical reality faced in the wireless 
multi-hop networks. Particularly, the interference caused by 
different transmitters may accumulate and is not binary, i.e., 
does not vanish beyond any specific border. Moreover, a 
node may successfully receive a message even when there 
are other concurrently transmitting nodes in its interference 
range. In [2], [3], it is argued that the performance of protocol 
model based algorithms is inferior to those with more realistic 
and fundamental physical interference model. More recently, 
Moscibroda et al. [4] show experimentally that the theoretical 
limits of any protocol, which obeys the laws of protocol inter- 
ference model, can be broken by a protocol explicitly defined 
for the physical interference model. However, although the 
physical interference model enjoy a much better performance, 
the scheduling problem under this model is notoriously hard 
to resolve. 

In this paper, we study the link scheduling problem under 
the physical interference model. First, we formulate the prob- 
lem as an Integer Linear Programming (ILP) problem, which 
is similar to the work in [5]. Since it is NP-complete [6], 
we get the approximate algorithm by using the randomized 
rounding method, and the probability bound of getting a 
guaranteed approximate ratio is given. Unlike the work in [5], 
we adjust the results of the rounding procedure based on all 
the constraints, and no iteration procedure is needed. 

Since the centtalized algorithm needs the global information 
of the network to do the scheduling, including the location of 
all nodes, the load of the links to be scheduled in a certain time 
slot, etc., collecting the required information may causes a lot 
of bandwidth overhead, especially when the network scale is 
large. Moreover, in some cases, the information is hard to get. 
The difficulties in the implementation of centralized algorithm 
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call for the distributed scheduling scheme that may achieve a 
fraction of the overall optimal performance. We then propose 
a distributed approach based on physical interference model. 
Each sending node implements the approach in three phases at 
each time slot, and all the constraints of the ILP problem in the 
centralized algorithm are satisfied. To analyze the performance 
of the algorithm, the approximation ratio is given. 

The rest of the paper is organized as follows. After giving 
an overview of related work in Section [HI the network and 
interference model is described in Section [III] In Section 
IIVI scheduling problem under physical interference model is 
formulated and we present the centralize algorithm to get the 
approximate solution. In section [V] a distributed approach is 
proposed. Section [VI] presents the simulation results of the 
proposed algorithmss and Section [VTIl concludes the paper. 

II. Related work 

The problem of scheduling link transmissions in a wireless 
network in order to optimize one or more of performance 
objectives (e.g. throughput, delay, fairness or energy) under 
interference constraints has been a subject of much interest 
over the past decades. 

As mentioned in the previous section, an accurate modeling 
of interference is fundamental to the scheduling problem. Most 
of the scheduling mechanisms proposed for wireless multi-hop 
networks use the protocol interference model, such as [7]-[l 1]. 
These algorithms usually employ an implicit or explicit color- 
ing strategy, and simply neglect the aggregated interference of 
nodes located farther away. However, the interference caused 
by a transmitter is not binary, i.e., it does not vanish beyond 
any specific border, and may accumulate amongst multiple 
concurrent transmissions. 

Only a few latest work have considered physical interfer- 
ence in the context [12]-[20], etc. In [12], Jain, et al. formulate 
the problem of scheduling under physical interference model 
as an LP problem. Unfortunately, no polynomial time solution 
and simulation-based evaluation of scheduling is given. [13] 
also provides an exponential-time LP formulation. In [14], 
Brar et al. present a heuristic scheduling method that is based 
on a greedy assignment of weighted links. Although it is 
based on physical interference model, the approximation factor 
of the algorithm is given only when nodes are uniformly 
distributed. The work of [15] considers physical interference 
in the minimum length scheduling problem. It uses a power- 
based interference graph model, which describes the interfer- 
ence relationship of every two links according to the SINR 
of the receiver. However, the model fails to consider the 
accumulation effect of interference. In [16], approximation 
algorithms for packet scheduling to minimize end-to-end delay 
with physical interference model are proposed. The works of 
[17], [18] study the problem of scheduling edges with SINR 
constraints to ensure that some property (e.g., connectivity) is 
satisfied. Similar with the work of [19], [20], they take power 
control, scheduling or routing together into account. 

The work of [5] is similar to our centralized algorithm 
discussed in section [IV] It also formulates the scheduling 



problem as ILP, but the randomized rounding procedure uses 
an iteration method which is more time-consuming than ours. 
Because the rounding procedure would not stop until a feasible 
solution is found, it costs a very long time to converge, and 
in some cases, it may not converge at all. 

Various distributed algorithms have been proposed for find- 
ing good approximations of the scheduling problem based on 
the protocol interference model (e.g. [21]-[23], etc.). Only 
a few previous work propose distributed algorithm based 
on the physical interference model. The work of [24]-[26] 
propose distributed algorithms based on the physical inter- 
ference model, which is lattice-throughput-optimal. But the 
approximation ratio is not given. The work of [27] develops a 
constant-time distributed random access algorithm for schedul- 
ing and gives the performance bound of the algorithm. In 
[28], a distributed and randomized protocol is proposed, which 
uses physical carrier sensing to reduce message overhead. It 
is similar to our approach discussed in section [V] But the 
protocol in [28] is sensitive to the scale of the network and 
spends much time to learn this information. 

III. Network and interference models 

We abstract a wireless multi-hop network as a directed graph 
G(V, E) where V is a set of vertices denoting the nodes 
comprising the network and £ is a set of directed edges 
between vertices representing inter-node wireless links. The 
Euclidean distance between any two nodes Vi,Vj G V, is 
denoted by d(vi,Vj), Let G E denotes the edge between 
Vi and Vj. 

Each node is assumed to be equipped with a single 
transceiver working in the half-duplex way, and all nodes share 
a common channel. So a node can not send and receive packets 
simultaneously. All antennas are omnidirectional. 

It is assumed that the network is using Time Division 
Multiple Access (TDMA) MAC protocol. The time is divided 
into slots of fixed length, which are grouped into frames. To 
increase capacity, spatial reuse TDMA (STDMA) [29] can 
be used, which is an extension of TDMA. The capacity is 
increased by spatial reuse of the time slots. 

We use the physical interference model [1] to describe 
the interferences between active links. In this model, the 
successful reception of a transmission depends on the received 
signal strength, the interference caused by nodes transmitting 
simultaneously, and the ambient noise level. The received 
power P r (si) of a signal transmitted by sender Sj at an 
intended receiver Vi is 

P r { Si ) = P(sj) • g(si,n), 

where P(si) is the transmission power of Si and g(si,ri) is 
the propagation attenuation (link gain) modeled as g(si, r^) = 
d(si,ri)~ a . The path-loss exponent a is a constant between 
2 and 6, whose exact value depends on external conditions of 
the medium (humidity, obstacles, etc.), as well as the exact 
sender-receiver distance. As common, we assume that a > 2 
[1]. Given a sender-receiver pair (sj,rj), we use the notation 
I r (sj) — P r (sj) for interference from any other sender Sj 
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concurrent to S{. The total interference I r at the receiver Vj 
is the sum of the interference power caused by all nodes 
that transmit simultaneously, except the intending sender S{. 
Mathematically, we have I r :— E s ev\{s ■} ^r(sj)- Finally, 
let N denote the ambient noise power level. Then, n receives 
transmission successfully from Si if and only if 

P r (8i) 



{v i ,v j )-o> x ij 



(1 - 4-)A 



Pjvk 



> 0, Ve y - G E,Vt (5) 



efcj6-B\{ei 3 -} d(v k ,Vj)- 



4 G {0,1}, Ve - G B,Vt 



(6) 



SINR{r,) 



N + J2 s . evx{Si} I r (sj) 

P(sj) ■ g(si,n) 

N + E S] ev\{ Sl } P ( s j) ■ 9(sj,n) 

>0, 



d(si ,ri) 



N 



where /3 is the minimum SINK threshold required for a 
successful message reception. 



The objective function of the formulation is to maximize the 
total network throughput, which is defined as the total traffic 
transmitted per slot. Constraint (Q]) guarantees that each active 
edge should be scheduled at least once. Constraint (O and 
constraint ([3]i make sure that each node can only receive or 
send signal from or to one another node. This is because each 
node has only one transceiver. Under constraint each node 
can not send and receive at the same time because it work in 
a half-duplex way. Constraint (O expresses the required SINR 
threshold that should be satisfied in order to have a successful 
reception at the receiver. The item (1 — x^)A ensures that the 
inequality is also satisfied when link is not scheduled in 
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In this section, we first formulate the link scheduling prob- 
lem under the physical interference model as an Integer Linear 
Programming problem, then give an approximate algorithm 
using randomized rounding method, which can be done in 
polynomial time. 

A. Problem formulation 

We denote by Vt C V and Vr C V the set of transmitting 
and receiving nodes respectively. The time is divided into slots 
of fixed length and a frame is composed of T times slots, the 
length of which is constant. In time slot t G [1, T], there exists 
at least one edge being scheduled to transmit. We denote by 
bij{ei,j G E) the total traffic rate through link and we 
use these to measure throughput. The goal of the scheduling 
method is to transmit all the edges in a frame to gain the 
maximum throughput. 

To formally formulate the problem, the boolean variables 
x\a is defined as 

1 if link eij is scheduled to be transmit in slot t 
otherwise 

Based on above assumptions, the linear integer program- 
ming formulation can be written as follows, 



1 T 

yE E bi 3 X ij 
t=l eyg-B 



s.t. J2 x ^ - 1; Ve ^ e E (1) 

t=l 

4<!> VvjeV R ,Vt (2) 

(3) 



J2 ^eV T ,Vt 
Vjev R 

E 4 + E 4 < v Vl g v s n v R , vt (4) 



B. Approximate algorithm with randomized rounding 

Since ILP problems are NP-complete, there is no efficient 
algorithm known for solving them in bounded time (and there 
can not exist any unless P = NP). The above formulation 
does not help us solve the scheduling problem. However, it 
does guide us to a natural relaxation which helps to find a 
good approximate algorithm. 

1) Randomized rounding procedure: We can use random- 
ized rounding method to find efficient and near optimal 
solutions [30]. In this procedure, each variable x\a will be 
relaxed to be in the range from zero to one and therefore 
converted from integer variables to fractional ones, i.e., xjj G 
[0,1], forVeij G E,Vt. The ILP problem is then relaxed 
to be a linear programming problem, which is called LP- 
relaxation. The optimal solution of the LP-relaxation can be 
treated as a probability vector with which we choose to include 
a link to a specific times lot. Denote the optimal solution as 
x\y In that sense, if is near one, it is likely that this 
link will be included in the current times lot. If, on the other 
hand, the optimal solution of the LP-relaxation problem is near 
zero, it will probably not be included in the current times lot. 
More specifically, the rounding procedure can be described as 
follows, 



1- 
0. 



with probability . 
with probability 1 — 



(7) 



However, after the rounding procedure, all the constraints 
may not be satisfied. For example, we suppose in constraint 
(O, for Vvj G V R ,Vt, there are two nodes v±,V2 G Vt- If 
= 0.3 and = 0.7, constraint (O is satisfied. But 



after the rounding procedure, x%j = 1 with probability 0.3, 
and X2j — 1 with probability 0.7. So J2 Vi ev T 4 = ^ w ^ 
probability 0.21. That is with probability 0.21, constraint (O is 
not satisfied. To overcome this problem, the work in [5] uses a 
iteration method. The rounding procedure would not stop until 
a feasible solution is found, which costs a very long time to 
converge. Even worse, it may not find any feasible solution 
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Algorithm 1 Centralized Approximate Algorithm 



Consider all links in E: 

Relax the integer linear programming problem to be a 

linear programming problem (LP-relaxation); 

Find the solutions of the relaxed LP {x\^ : Ve^ G E,t G 

T}; 

Do randomized rounding and get the rounding solution 
{x lj : Vey e E,t eT} from equation Q; 
for each t G T do 



Get the set 0* = {c 
Consider all links e 
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= i}; 

( in non-decreasing order of 



if constraint [2] or constraint [3] or constraint [4] is not 
satisfied then 

4- - 0; 

Move eij away from </>'; 
end if 

Consider all links G in non-decreasing order of 

if constraint [5] is not satisfied then 

x% = 0; 

Move away from </>*; 
end if 
end for 

Consider all links ey G -E: 
for each link en E E do 



4=1 ^ij 



< 1 then 



Consider the biggest element i' ? for all t G T: 

t 



1; 



Execute from line (5) to (17) 
end if 
end for 



in the end. Unlike the work in [5], we adjust the results 
of the rounding procedure based on all the constraints, and 
no iteration procedure is needed. The detail of our proposed 
algorithm is discussed in the following part. 

2) Approximate algorithm: we denote by $ — {e.y : x\^ — 
1} the set of edges that will be transmitted after the rounding 
procedure at time slot t. As discussed above, all the constraints 
of the ILP program may not be satisfied. So the solution of the 
rounding procedure will be adjusted in algorithm 1 to satisfy 
all the constraints. As the solution of the relaxed LP problem 
{x\j : Vey G E, t G T} is the probability for the wireless 
links to transmit, we reorder the elements in according to 
the probabilities in the non-decreasing order. 

Algorithm [T] proceeds in two phases: phase 1 is correspond- 
ing to line 1 to 4; while phase 2 is corresponding to the rest of 
the algorithm. In phase 1, the ILP problem is relaxed to be a 
linear programming problem which is solved by the random- 
ized rounding procedure. As mentioned before, there exists a 
problem that the rounding solution \x\j : Vey G E, t G T} 
may not satisfy all the constraints in the ILP problem. It can 
be resolved in phase 2. From line 5 to 17, it is checked slot by 



slot whether constraints (2) to (5) are satisfied. The solution 
of the relaxed LP {x\a : Vey G E, t G T} can be viewed 
as the probability of the transmission of the links. The links 
with larger value of x\^ have higher priority to transmit. So 
we check the links in a non-decreasing order of x\y In each 
slot, if a link doesn't satisfy any of the constraints of (2), (3) 
or (4), it will not transmit in the current slot. After all links 
in are checked, constraint (5) is checked from line 12 
to 16. Because constraint (5) is the interference constraint, 
with which the transmission of a link is related to all the 
other transmitting links, it is checked after constraint (2) to 
(4) being checked. Constraint (1) guarantees that each active 
link should be scheduled at least once in a frame. However, 
after the procedure of randomized rounding and the check of 
constraint (2) to (4), x\: of a link may be rounded to 0, or 
adjusted from 1 to 0, leading that constraint (1) is not satisfied. 
From line 18 to 25, constraint (1) is check. For a link that 
doesn't satisfy constraint (1), find the time slot t in which 
it has the largest probability to transmit (line 21). Then let 
Cij has probability 1 to transmit, which guarantees that e y - is 
scheduled once. Constraint (1) is then satisfied. 

3) Performance analysis: In this part, the performance of 
algorithm ?? is analyzed, including the complexity and the 
approximation ratio. It can be seen that algorithm Q] can finish 
in polynomial time, and the probability lower bound of being 
a (1 — (^-approximation algorithm (V0 < 8 < 1) is calculated. 

The complexity of phase 1 depends on the algorithm solving 
the LP-relaxed problem. Let it be denoted by 0(P). Since it 
is a linear programming problem, it runs in polynomial time. 
Thus P is polynomial. Let TV be the total number of the links 
in E. It can be easily seen that the complexity from line 5 to 17 
in algorithm 1 is O(NT), because there are two loops there. 
Similarly, the complexity from line 18 to 25 is 0(N 2 T). Thus 
the complexity of the whole algorithm is 0(max(P, N 2 T)), 
which is polynomial. 

The approximate ratio of algorithm ?? can be got from the 
following theorem: 

Theorem 1 (Approximate ratio of the centralized algorithm): 
For V0 < 6 < 1 and -6 < AA/A < 1 - 6, the probability of 
algorithm ?? being (1 — #)-approximate to the optimization 

is lower bounded by 1 — e 2 , where A is throughput 
calculated by the LP-relaxation, AA is the variation of the 
throughput in phase 2. 

Proof: In phase 1, let A ran( i and A be total throughput 
calculated by the randomized rounding in phase 1 and LP- 
relaxation respectively, i.e. A rand = Y2~2t=i2~2 ei eE b ij x ij' 
and A = ?f2~2t=i 2~2 eij eE b ij^ij- F rom equation (0, we can 
get E{x\ j ) = x\y Thus, 



1 T 

E(A rand ) = E(-J2Y1 b V x W 

t=l e i:j eE 



t=l ea&E 



5 



.4 



(8) 



According to equation ^ and Chernoff Bound [30], we can 
get the following bound, 

Pr(A rand > (1 - S)A) = 1 - Pr(A rand < (1 - 5) A) 

> 1 - < 5 < 1 (9) 

Let A be the final solution of algorithm ?? and Ao P t be the 
optimal solution of the ILP problem. Obviously, 



A < A 0pt < A 



(10) 



In phase 2, the randomized rounding solution is checked 
whether it satisfies all the constraints in the ILP problem, 
and some adjustments are made, leading that the throughput 
changes. Suppose the throughput changes by AA after phase 
2, and it follows that 



I i — . 1 , 



AA 



(11) 



Substituting equation (fTTb into (O, we can derive 

Pr(A> > 0<6< 1 (12) 

A 

In phase 2, the worst case is that all links are adjusted not to 
transmit, leading the throughput A = 0, i.e. A ran d + AA = 0; 
while the best case is that after the adjustment, the solution 
approaches to the optimal, leading A = A ran d + AA < A. 
Thus we can get the variation range of AA as follows, 



Arand < AA < A - A 



rajid 



(13) 



Let 6 = S — ^4-. Since < S < 1, it can be derived that 

A 

—8 < AA/ A < 1 — 8. Equation (TT2l) can be written as follows, 

Pr(A > (1 - 8) A) > 1 - e 5 , (14) 

where -1 + < 8 < 1 + (from equation (O), and 

-8 < AA/A<\-8. 

Because (0,1) C [-1 + Ara £ d , 1 + A ™£ d ], according to 
equation ( fTotlfPfl) . we can get 

Pr(A > (1 - 8)A 0p t) > Pr(A > (1 - 6)A) 

>l-e * (15) 

where 8 G (0, 1), and -8 < AA/ A <l-8. 



V. The distributed algorithm 

The approximate algorithm discussed in section [TV] uses the 
physical interference model to represent the interference and 
has a good performance. However, as mentioned in section 
[U the centralized algorithm requires the global information 
of the network to do the scheduling. Collecting the required 
information may causes a lot of bandwidth overhead, espe- 
cially when the network scale is large. Moreover, in some 
cases, the information is hard to get. The difficulties in the 
implementation of centralized algorithm call for the distributed 



scheduling scheme that may achieve a fraction of the overall 
optimal performance. 

In this part, we extend the centralized algorithm ?? to 
a distributed solution, which takes the physical interference 
model into account and implements in three phases. It is 
assumed that nodes in the network can perform physical carrier 
sensing, and they can set the carrier sensing range to different 
values. As we show in the following part, by properly tuning 
the carrier sensing range, all the constraints of the ILP problem 
in section [IV] can be satisfied. The approximation ratio is also 
proved. 

A. Carrier sensing range calculation 

Let d m ax and d m i n be the maximum and minimum link 
length in the network respectively. The length diversity k is 
defined as 

k = [log (d max / dmin) \ (16) 

Suppose each node transmits at the same power P. As the 
ambient noise power is much less than the interference power, 
it is omitted here. The following theorem describes how to tune 
the carrier sensing range for each node. 



\ 



\ 



/ 



VI/ Vj \ 




Fig. 1. The Euclidean plane is divided into a series of rings 

Theorem 2 (Carrier sensing range): If the carrier sensing 
range is set to be Rc = p2 k , the transmitting links in the 
network will not be interfered by all the other links, where p 
is defined as 

p = 4(2Tr[3^±)i, (17) 
a — 2 

a is the path-loss exponent and (3 is the minimum SINR 
threshold). 

Proof: Without loss of generality, suppose link G E 
is transmitting data from node Vi to Vj. we normalize the 
minimum link length d m in to be 1 . From the definition of the 
length diversity, we can get 2 k < d ma x < 2 k+1 . The length 
of link eij satisfies < d ma x < 2 k+1 , thus the perceived 
power at Vj from Vi is at least 



P 

p. > 

■M) - 2 (*+i) 



(18) 



The plane can be divided into a series of rings with the center 
located at sending node V{, which can be seen in figure Q] 
They are denoted as ring 0, 1, • • •, respectively. Since the 
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region of ring is in the carrier sensing range of the sending 
node Vi, there are no other sending nodes locating in ring 0. 
For a sending nodes v m located in ring 1, it must be at least 
d m j > Rc — 2 k+1 away from Vj. Since each sending node 
must be outside the carrier sensing range of other sending 
nodes, there are at most 7r/arcsin(5) = 6 sending nodes can 
transmit concurrently in ring 1. Consequently, the aggregated 
interference at is 



E^ 

i=i 



A p 



1=1 a io 



< 



6P 



(R c -2 k + 1 ) a 

In ring 2, the sending nodes are at least 2Rc — 2 k+1 away from 
Vj. there are at most 7r/arcsin(i) < jjj = 4tt sending nodes 
can transmit concurrently. So the aggregated interference at Vj 
is 

E p ^ 



{2R C - 2 k+1 ) a 



Similarly, in ring to, the sending nodes are at least mRc — 
2 k+1 away from Vj. there are at most tt/ arcsin( < 
1 Z m — 2to,7t sending nodes can transmit concurrently. So 
the aggregated interference at Vj is 

2m7rP 



2m7T 

E 

1=1 



Pij< 



(mR c - 2 k + 1 Y 



Consequently, the accumulated interference at the receiving 
node vj is upper bounded by 



I„ 



< 



E 



2mnP 



< x (mRc - 2 k+1 Y 
2mirP 



2 *+n 



< 



E ( mp 2 k 

m=l v 
27rP y, 

2ttP to 
2 ka 2-, (mp/2) c 



< x (mp - 2Y 



(19) 



2ttP 



.a ^— < 



< 



2' £Q ( / o/2) Q " to- 
2ttP a - ] 



2 (fe-l)a p a Q _ 2 ( 20 > 

where (fT9l follows because x — 2 > x/2, Va; > 4 and p > 4, 
given that /3 > 1 and a > 2; and (f20b follows from a bound 
on Riemann's zeta function [6]. From equation ( TPTI i. (fl~8T > and 
( |20] >, we can get the SINR at the receiving node Uj as follows: 

s/ivp,,, = -ii- 



> 



2( fc +i)° 



2ttP 



q-1 



2 (fc-i)apQ Q -2 

/3 



B. Algorithm description 

The frame structure is shown in Figure [2] A frame is 
composed of fixed length time slots. A time slot is divided into 
3 phases, namely, carrier sensing phase, RTS-CTS phase and 
data-ack phase. Our distributed scheduling approach allows 
each node to implement it in each of the three phases. Next, 
we will present the details of the approach and show that after 
the three phases, all the constraints in the ILP problem are 
satisfied. 

Frame 



Slotl 


Slot 2 


Slot 3 












^ 




Carrier 










RTS 


CTS 


Data 


Ack 


sensing 











Fig. 2. Frame structure in distributed approach 

1 ) Carrier sensing phase: The purpose of this phase is to 
guarantee that there are no other nodes sending data in the 
carrier sensing range of a transmitting node. Each node in 
the network sets its carrier sensing range to be Rc, which is 
calculated in theorem [2] So according to theorem [2] the node 
will transmit its data without any interference. Constraint (0 
in the ILP is satisfied. 

Consider a transmission slot starts at time t. The carrier 
sensing phase begins at the beginning of each time slot. It lasts 
for a period of t c . Each node that wants to transmitting data 
randomly selects a time t s and sends a SENSING signal, which 
lasts for a period of t (i < t s < t + t c — r). The purpose of the 
SENSING signal is to let each node sense whether there are 
other transmitting nodes in its carrier sensing range. If a node 
do not sense any signal before sending the SENSING signal, 
it will occupy this time slot, and enter phase 2; otherwise, the 
node randomly selects another time slot and waits to transmit 
again. So the nodes located in each other's carrier sensing 
range will transmit in different slots, and no interference will 
occur. For each node, the purpose of randomly selecting a time 
to send SENSING signal is to prevent that all nodes send the 
SENSING signal concurrently, and they all will sense nothing. 

If a link is scheduled, the sending node will not send the 
SENSING signal again until it can sense nothing, i.e., all the 
sending nodes that have been scheduled at least once will wait 
for the nodes that haven't been scheduled to transmit. This rule 
guarantees that all the links are scheduled at least once in a 
frame. So constraint (fl~|i is satisfied. 

2) RTS-CTS phase: In this phase, the sending node first 
sends RTS signal to the corresponding receiving node. If the 
receiving node receives RTS, it will respond CTS signal to the 
sending node. If the sending node receives CTS signal, it will 
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enter phase 3. Otherwise, it will randomly select another time 
slot to transmit. 

This phase guarantees that constraint (O in the ILP is 
satisfied, i.e., a node can not receive data from more than 
one nodes simultaneously. An example can be seen in figure 
[3] Node A and B are both transmitting to node C. Since node 
B is outside the carrier sensing range of node A, A and B 
can transmit simultaneously in phase 1 without interference. 
Because C can not receive data from more than one node 
simultaneously (constraint |2j, the data from A or B will not 
be received. Through the RTS-CTS procedure, A or B will 
transmit in another slot, and constraint (f2]l is satisfied. 



v 



defined in theorem |2] and d max is the normalized maximum 
link length). 

Proof: It can be seen that each sending node in the 
network can be located in a certain circle with radius being 
the carrier sensing range Rc- Without loss of generality, we 
consider a circular area Q Vi with a sending node v.- L located in 
the center. Node vi is sending data to vj. Assume an optimal 
algorithm OPT can schedule at most q links in Topt = 1 slot 
whose sending node is located in $ Vi . According to theorem 
|2] to scheduling these q links, our proposed approach needs 
T < q time slots. So the approximation ratio is 



T 



T, 



OPT 



<q 



Next, we will calculate the maximum value of q. The best 
case is that there are no node sending data outside $ Vi , and 
all the links in $ Vi can transmit successfully. For node Vj, the 
perceived SINR level is 



SINR{vj) 



< 



Fig. 3. An example in phase 2 

3) Data-ack phase: In this phase, the sending node sends 
data to the receiving node, and waits for the ack signal from the 
receiving node. If the receiving node receives data correctly, it 
sends back ack signal to confirm the correct reception. If the 
sending node doesn't receive the ack signal, it will transmit 
the data again in the next slot. 

Because of the distribution feature of our approach, a node 
can decide not to send data to more than one receiving nodes 
simultaneously. So constraint (01 is satisfied. Usually, the 
carrier sensing range is larger than the transmission range. 
Thus, when a node is transmitting, no other nodes will transmit 
to it. In addition, because of the half duplex character of 
wireless nodes, a node will not send data to others when it 
is receiving data. So constraint <j4j is satisfied. In phase 1, 
constraint ([TJ and (O are satisfied. Constraint (0 is satisfied in 
phase 2. So all the constraints are satisfied in our distributed 
approach, which can be viewed as a distributed solution of 
the ILP. Next, we will analyze the approximation ratio of the 
distributed approach. 

C. Performance analysis 

The approximation ratio of the distributed approach can be 
described in the following theorem: 

Theorem 3 (Approximation ratio of the distributed approach) 
The approximation ratio of the distributed approach is at 
most d" lax (p + 2) a /(3, i.e., the distributed approach is an 
0(d" lax ) approximate algorithm (p, a and (3 are the same as 



p 

qP/(p2 k + 2 k + 1 ) a 
= SINR M ax(vj) 

where (fSTJ follows because dij > 1 and dkj 
According to SINRmax{vj) < fl, we can get 

q < 2 ka (p + 2) a /0 
< C,(P+2)7/5 
= 0{d a max ) 



(21) 



< 2 fe+1 . 



(22) 



where (l22l follows from the definition of diversity dl6) . ■ 
From theorem [3] we can see that the performance of the 
distributed algorithm is closely related to the link length 
diversity. The approach performs well when the diversity is 
small. It can be seen distinctly from simulation (section I VP , 

VI. Simulations 

In this section, we present the simulation results that verify 
the performance of our proposed approximate algorithms in 
this paper. 

A. Simulation scenario 

We create a scenario where wireless nodes are uniformly 
distributed in a 100 x 100 square. In the simulation, we 
normalize the transmission range of a node to be Rt = 1, 
and the interference range to be Ri = 2.5. For simplicity, all 
nodes are supposed to transmit at the same power level. We 
also assume N = —90dBm, (3 — lOdB, and a = 4, which are 
similar with [2] and [24]. Let the frame length be T = 100. 
.The whole simulation runs 100 times. In each simulation, n 
communication pairs are randomly selected. Suppose each of 
the selected links has the same traffic rate which is normalized 
to be 1. 
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Number of communication pairs n Approximation ratio (1-6) 



(a) (b) 




Approximation ratio (1-0) Approximation ratio (1 -6) 

(c) (d) 
Fig. 4. The performance of the centralized algorithm 

B. Centralized algorithm 

We first compare the throughput performance between our 
proposed algorithm and several other centralized algorithms. 
The algorithm in [10] is based on protocol interference model 
and a polynomial-time coloring algorithm is used. To express 
conveniently, we call it Algorithm Protocol Model (PM). In 
[14], the algorithm is based on physical interference model, 
which is the same as ours. It is a heuristic scheduling method 
that is based on a greedy assignment of weighted links. We 
call it Algorithm Physical Greedy (PG). The authors in [15] 
construct the conflict graph using the physical interference 
model, which does not consider the accumulation effect either. 
It is called Algorithm Physical Conflict Graph (PCG). The 
throughput comparison of the several algorithms is shown in 
Fig |4(a)| The notation OPT Bound here stands for the optimal 
solution of the LP-relaxation. The notation OPT is the optimal 
solution of ILP, which is got by traversal searching. Algorithm 
APP implies our approximate algorithm. 

It can be seen that our proposed APP algorithm is well 
approximate to the optimal solution of ILP, which is a little 
lower than the upper bound. It also outperforms the protocol 
interference model based algorithm (PM algorithm) and other 
physical interference model based heuristic algorithms (PG 
and PCG algorithms). It performs ever better when more 
nodes begin to transmit. The PM algorithm performs the worst, 
because it uses the protocol interference model, which is less 
precise than the physical interference model. When more and 
more nodes begin to transmit, the throughput drops for all the 
algorithms. This is because when interference is more severe, 
more transmissions begin to conflict with each other. 

Then we evaluate the influence of the parameter A A/ A 
(defined in theoremHJ to the (1 — #)-approximate lower bound. 
Here we denote A A/ A by A A/ A in the legends of Fig. 4. 
It can be seen from Fig |4(b)| that the lower bound decreases 



while A A/ A decreases, which is from the fact that more 
adjustments are made in the phase 2 of the algorithm Q] In 
our simulation, AA/A is rarely smaller than -0.3, implying 
that the approximation ratio is larger than 0.5 with probability 
larger than 0.7, and the approximation ratio larger than 0.4 
with probability larger than 0.9. 

When more links begin to transmit, more adjustments are 
made in algorithm [U i.e., A A decrease. The influence of the 
network scale n to the performance can be seen in Fig |4(c)| and 
|4(d)| We can get the following conclusions when n increases: 

1) If AA/A doesn't decrease, the probability lower bound 
increases. In Fig |4(c)| AA/A is fixed to be 0. 1 . The probability 
lower bound increases while more links begins to transmit. If 
AA/A increases, the probability lower bound also increases 
(see plot (n = 60, AA/A = -0.8)and (n = 90, AA/A = 0) 
in Figggj). 

2) If AA/A decrease, the probability lower bound dereases. 
It can be seen in Fig |4(c)| i. 

C. Distributed algorithm 

In this simulation, nodes are not uniformly distributed, 
which is different from the scenario in the centralized algo- 
rithm. We select the transmission pairs properly to control the 
link length diversity k. 

In Fig. [5] the throughput of the centralized algorithm and 
distributed approach are compared. It can be seen that the 
distributed algorithm approaches the centralized algorithm 
well when k = 0. The throughput of k — and that of k = 1 
are very close to each other; while it decreases sharply when 
k = 2. It is because the approximation ratio is 0(2 ka ), which 
is exponential to k. However, in real scenario k is not very 
large, typically k < 2. The distributed algorithm proposed in 
the paper is preferred to be implemented in the scenario that 
the link length diversity is small. The ideal scenario is that the 
nodes are located uniformly. 




Number of conmmunication pairs n 
Fig. 5. The performance of distributed algorithm 

VII. Conclusion 

In this paper we study the link scheduling problem with 
the physical interference model, which is more accurate than 
the protocol interference model but notoriously hard to handle 
with. The scheduling problem is formulated to be a ILP 
problem, and we propose our centralized algorithm based 
on randomized rounding. We analyze the performance of 
the algorithm and give the probability bound of getting the 
guaranteed approximation ratio. As a distributed solution of 
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the ILP problem, we also propose a distributed approach, 
which uses physical carrier sensing and implement in phases. 
The approximation ratio is also given, which is associated with 
the link length diversity. 
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